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We performed inelastic neutron scattering experiment to study magnetic excitation of O2 
molecules adsorbed in microporous compound. The dispersionless excitation with characteristic 
intensity modulation is observed at huj — 7.8 meV at low temperature. The neutron cross section 
is explained by spin dimer model with intradimer distance of 3.1 A. Anomalous behaviour in 
the temperature dependence is discussed in the context of enhanced magnetoelasticity in the 
soft framework of O2 molecule. 
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Oxygen molecule, the second most abundant con- 
stituent in air, is essential to human beings. One of 
the interesting physical properties is that the diatomic 
molecule has magnetism with effective spin of S = 1. 
In solid state at low temperature the oxygen behaves as 
insulating low-dimensional antiferromagnets 1 thorough 
exchange interaction induced by the significant overlap 
of the molecular orbital. Under pressure the change of 
electronic structure brings exotic phenomena such as 
metallization, 2 superconductivity 3 and exotic magnetic 
orders. 4 ' 5 Now the understanding of magnetic oxygen 
as an elementary substance has been progressed. The 
next challenge is to design a magnet by treating oxygen 
molecule as a magnetic entity. Enterprising trial for the 
O2 based magnet was performed in adsorbing O2 mono- 
layer on graphite substrate. 6, 7 Neutron diffraction and 
magnetic susceptibility study revealed the structural and 
magnetic variation of the artificially designed triangular 
lattice of O2 molecule with the change of the oxygen den- 
sity and temperature. However the amount of the oxygen 
in the monolayer is too small to be detected for a dynam- 
ical probe such as neutron inelastic scattering technique. 
For further magnetic study availability of bulk size mag- 
net has been required. 

One solution is to utilize microporous compounds 
in a recently developing chemistry area, coordina- 
tion polymers science. 8-10 Particularly interesting is 
a copper coordination polymer {[Cu2(pyrazine-2,3- 
dicarboxylate)2(pyz)]2H20}„ abbreviated as CPL-1. 11 
Molecule formula is Ci6Hi2Cu2N60io- The space group 
is monoclinic P21/c and the lattice constants are a = 
4.693 A, b = 19.849 A, c = 11.096 A, and = 96.90°. 
The compound includes ID channel for physical adsorp- 
tion of molecules with the dimension of 4.0 A by 6.0 A 
in crystallographic a direction (Fig. 1(a)). A strong con- 
finement effect due to the restricted bore size enables the 
formation of ladder- like array of adsorbed oxygen 12 thor- 
ough van der Waals interaction. The intermolecular dis- 
tance in the rung/dimer direction is 3.28(4) A. Unlike the 
surface adsorption on graphite O2 molecules penetrate 
inside and the amount of the adsorption is one molecule 
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Fig. 1. (Color online) Crystal structure of the microporous cop- 
per coordination polymer (CPL-1). (a) The framework of non- 
adsorbed CPL-1. ID channel for molecules adsorption runs in 
crystallographic a direction, (b) Structure of CPL-1 with O2 ad- 
sorption. 



per one copper atom. Hence the O2 based magnet of bulk 
size is realized. The magnetic susceptibility and magne- 
tization measurements suggest a non-magnetic ground 
state with a spin gap. 12,13 The magnitude of the gap 
was roughly estimated as 88 K ~ 100 K (7.6 meV ~ 
8.6 meV). In this letter we report the first observation 
of the magnetic excitation in the prototype of artificial 
O2 magnet realized in CPL-1 by using inelastic neutron 
scattering technique. 



2 J. Phys. Soc. Jpn. 



Letter. 



Author Name 



Neutron scattering experiment was performed using 
PONTA spectrometer installed in 5G beamline in JRR- 
3M reactor in Japan Atomic Energy Agency. A series 
of scans were collected by changing incident neutron en- 
ergy for fixed monitor counts in a low-efficiency detec- 
tor between pyrolytic graphite (PG (002)) monochroma- 
tor and sample. 80' sollar collimator was used between 
monochromator and sample. Radial collimator was in- 
stalled between sample and PG horizontally focused an- 
alyzer. The energy of the final beam is fixed at f 4.7 meV. 
PG filter was placed after sample to eliminate higher 
harmonics. Projected full-width at half-maximum energy 
resolution of this configuration at elastic position is 0.96 
meV. 

Schematic view of O2 induction apparatus is shown in 
Fig. 2. ORANGE type cryostat is used to achieve T — 
2.0 K. 12 g of powder sample was prepared by the pre- 
viously reported method. 11 The sample is packed in Al 
cylindrical can (20 mm in diameter and 67 mm in length) 
which can be attached to 2 induction probe designed 
for the cryostat. VCR fitting with filter gasket (Swagelok 
Co. Ltd.) is used for the connection to stainless made O2 
induction tube. Spacer with the height of 1cm is put in 
the bottom of the Al can for the reservoir of extra bulk 
02- The adsorbed water molecules in the microporous 
channel are eliminated by heating the sample up to 400 
K in reduced pressure. 

Since it is difficult to synthesize the deuterated CPL-1 
we use protonated sample. To measure background from 
the host CPL-1 including large incoherent scattering of 
hydrogen atoms we performed precedent scans in the 
non-adsorbed CPL-1 for all scans. Afterwards the mea- 
surements for the O2 adsorbed CPL-1 were performed. 
O2 adsorption is processed by warming up the sample 
and by keeping the temperature at 95 K for an hour 
with O2 pressure of 80 kPa. 

Typical constant wave number (Q) scans of O2 ad- 
sorbed and non-adsorbed CPL-1 arc shown in Fig. 3 (a). 
The contribution of the adsorbed oxygen is detected as 
the enhanced intensity at the energy transfer (fko) of 
about 7.8, 10.6 and 13.7 meV. We assume that the inten- 
sity in O2 adsorbed sample is, as zeroth approximation, 
simple sum of the excitation from the adsorbed O2 and 
the background from the host compound. Then the sub- 
tracted excitation of the adsorbed oxygen is shown in 
Fig. 3 (b). Solid curves are the fit to the data by triple 
Gaussian functions. Among the observed three peaks the 
peak at fko = 7.8 meV shows typical behaviours of the 
excitation of a magnetic cluster: the peak energy does 
not change with Q and the peak intensity decreases with 
the increase of Q. In addition the peak width is within 
resolution limit as is indicated by shaded area. On the 
other hand the intensity of the peaks at fou) = 10.6 and 
13.7 meV is almost Q independent. These peaks can be 
ruled out from the magnetic origin since a magnetic neu- 
tron cross section is reduced at high Q due to magnetic 
form factor. 14 For some unknown reasons the subtracted 
intensity has negative values at 6meV < hu; < lOmcV 
and we presumed negative background for the Gaussian 
fittings. Probably O2 adsorption suppresses the incoher- 
ent scattering of hydrogen in host compound. The en- 
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Fig. 2. Schematic view of O2 induction apparatus. 



ergy scans after the negative background subtraction are 
shown in Fig. 3 (c). The characteristic behaviour of the 
magnetic cluster is observed more clearly. Q dependence 
of nuetron cross section at fku = 7.8 meV is shown in 
Fig. 4. Filled circles are the net magnetic excitation of ad- 
sorbed O2 obtained by similar data analysis in Fig. 3(c). 
Open circles are the scaled intensity of constant Q scans 
obtained in Fig. 3(c). 

The most probable spin model for the magnetic cluster 
is Hcisenberg S = 1 spin dimers from crystallographic 
consideration. 12 The effective spin Hamiltonian is de- 
scribed by H = —2JS\-S2- The calculated eigenenergies 
are E(S) = -J (S (S + 1) - 3/2) where S is the magni- 
tude of the total spin operator S — Si + S2 (5 = 0, 1, 
2). Each energy level is 25*+ 1 - fold degenerated and the 
energy levels satisfy Lande-interval rule. The magnetic 
excitation at low temperature is dominated by the tran- 
sition between non-magnetic singlet ground state and S 
= 1 triplet state. The neutron cross section for the mag- 
netic scattering is obtained 15 
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Here p\ (T) is a temperature factor that is proportional to 
the thermal population of singlet ground state: Pi(T) = 
1/ (l + 3aip(2J/k B T) + 5atp{6J/k B T)). F(Q) is mag- 
netic form factor of oxygen molecule, Q is a scatter- 
ing wave number and R is intradimer distance between 
oxygen molecules. For the calculation of F(Q) the den- 
sity of atomic P± orbitals is assumed for the spin 
density of an O2 molecule. 16 Then it is obtained that 
F (Q) = Jo°° dr So d6ri sin3 0e ~ 2bT ' 2 [cosh(26ri? cos 6) - 
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Fig. 3. Excitations of O2 based magnet from inelastic neutron 
scattering measurements at T = 2.0 K. (a) Constant Q scans at 
Q = 1.5 A -1 in O2 adsorbed CPL-1 (filled circles) and non- 
adsorbed CPL-1 (open circles), (b) The subtracted intensity. 
Shaded area is the experimental resolution, (c) Net magnetic 
excitations in narrow energy range at various Q's (offset for clar- 
ity). 

1 ] e »Qr cos(0-/3) Here Ro = ! 21 A and b = 4.1 A~ 2 . 17 R 
is the internuclear distance, b is a fitting parameter, and 
[3 is the angle between molecular axis and scattering vec- 
tor Q. F(Q) is numerically obtained by powder average 
F(Q) = fF(Q)dQ. 

The Q independent magnetic excitation at 7.8 meV in 
Fig. 3(c) is consistent with eq. (1) with 2J = —7.8 meV. 
Characteristic feature of the dimers neutron cross section 
is identified in constant energy scan at hw = 7.8 meV in 
Fig. 4. Enhanced peak structure at Q ~ 1.3 A^ 1 and 
drastic suppression of the intensity with the increase of 
Q are perfectly reproduced by the calculation based on 
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Fig. 4. Wave number dependence of the neutron cross section of 
the adsorbed magnetic O2. Symbols are the subtracted data and 
solid curve is the calculation of S = 1 spin dimcr model in eq. 
(1) in the text. 

5=1 spin dimer model with R = 3.1 A. The value of R 
is reasonably consistent with intermolecular distance of 
R = 3.28 A obtained by synchrotron x-ray diffraction. 12 
These facts mean that magnetic excitation at T = 2.0 
K is explained by the adsorbed O2 dimers. Furthermore, 
the energy scale is consistent with the spin gap estimated 
by bulk measurements 7.6 ~ 8.6 meV. 12 ' 13 

Constant Q scans at various temperatures is shown 
in Fig. 5(a). The magnetic peak at fuo — 7.8 meV is 
suppressed with the temperature and almost disappears 
a t T > 50 K. As is shown in Fig. 5(b) the peak inten- 
sity goes down more rapidly than the temperature factor 
pi(T) shown by the solid curve. The data at T < 40 K 
is rather consistent with two level system with singlet 
ground and octet excited states described by an empiri- 
cal formula p[(T) = 1/ (1 + 8 exp(2J/fc B T)). 

Since the framework of O2 molecules formed by van 
der Waals interaction is flexible, a magnetoelastic effect 
can be enhanced enough to affect the energy spectrum 
of the effective spin model. The calculation of the inter- 
molecular potential of O2 dimer showed that the energy 
level, particularly in high energy range, depends both on 
spin state and geometrical configuration. 18 The ground 
state is S — singlet with H configuration in Fig. 5 (c) 
and the first excited state is S = 1 triplet with the same 
H configuration. However, the second excited state is S 
= 2 quintet with X configuration and the energy interval 
from the first excited state is roughly — 2J. This value 
is two times smaller than that by Lande-interval rule of 
the spin dimer model. These calculations indicate that 
the effective spin Hamiltonian works only for the low 
energy excitation and the inclusion of van der Waals po- 
tential is essential to the higher energy excitations. For 
the understanding in quantitative level the consideration 
of specified potential including the interaction between 
O2 and CPL-1 is required in the calculation. 

In summary we report the magnetic excitation of ar- 
tificial O2 magnet in microporous compound CPL-1 by 
using inelastic neutron scattering technique. While the 
neutron cross section at low temperature is completely 
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Fig. 5. Temperature dependence of the magnetic excitation, (a) 
Constant wave-vector scan at Q = 1.5 A -1 at various tempera- 
tures (offset for clarity). Solid curves are guide for the eyes, (b) 
Temperature dependence of the integrated intensity. Solid and 
dashed curves are temperature factors described in the text, (c) 
O2 dimer in H configuration (left) and X configuration (right). 



explained by existent spin model, the temperature de- 
pendence shows anomalous behaviour due to enhanced 
magnetoclastic effect of the adsorbed O2 molecules. Our 
study demonstrates that the conventional dynamical 
probe for physics study, inelastic neutron scattering tech- 
nique, is powerful tool for an O2 based magnet in recent 
developing chemistry area. 
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